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Abstract 


A previously developed technique allows an estimate of integral 
mixing to be obtained from an image of laser scattered light from par- 
ticle seeded fuel in the hypervelocity flow through a scramjet combus- 
tor. This previous mixing analysis formulation contains an assumption 
of a constant velocity flowfield across the plane of the fuel plume 
image. For high-speed scramjet combustors, the velocity flowfield is 
quite uniform and an assumption of constant velocity works well. 
Applying this same mixing analysis technique to fuel plume images 
obtained from a mid-speed scramjet combustor makes it desirable to 
remove the constant velocity assumption. This is due to the non-uni- 
form velocity flowfielcls present in mid-speed scramjet combustors. A 
new formulation of the mixing analysis methodology is developed and 
presented so that the technique can be applied to a mid-speed scramjet 
combustor without the need to assume a constant velocity flowfield. 

Nomenclature 


A Area 

f Fuel-air mass ratio 

H Fuel injector ramp height 

I Image intensity 

m Mass flow rate 

N Number density 

N a Avogadro’s number 

U Local velocity 

Subscripts 


W Molecular weight 

X Distance from fuel injection location 

Y Mass fraction 

ij M Mixing efficiency parameter 

p Mass density 

t Image exposure time 

<j) Fuel equivalence ratio 


c 

Camera 

ROI 

Region of Interest 

f 

Fuel 

RX 

Reactable 

L 

Laser 

ST 

Stoichiometric 


Introduction 

The objective of scramjet combustor ground tests is to acquire data useful in quantifying tire perfor- 
mance of the modeled scramjet. Key to performance is the fuel mixing achieved within the combustor, 
which can be deduced from flowfield data. One way to achieve this result is to acquire fuel concentra- 
tion distributions within the fuel plume stmcture through the use of fuel plume imaging. Image data 
obtained in this manner can be processed to yield an estimate of the integral fuel mixing. Initially this 
technique was applied to high-speed scramjet combustors and tire analysis method utilized an assump- 
tion of constant velocity across the plane of the fuel plume image [1], The purpose of the current work is 
to extend the range of applicability of fuel plume image mixing analysis to include mid-speed scramjet 
combustors, for which the velocity flowfield is much less uniform and cannot be assumed to be con- 



stant. A reformulated fuel plume image mixing analysis method is developed and presented along with 
guidelines to help determine which method of analysis is appropriate. 


Fuel Plume Imaging 


Fuel plume imaging relies on Mie scattering from particle seeded fuel which is illuminated by a 
sheet of light optically formed from a laser beam. References 1 , 2, and 3 contain specific details con- 
cerning optimization of operational parameters for the purpose of obtaining clean fuel plume images. 
Basically, silica spheres are used to seed the fuel prior to injection. It is assumed that the seed particles 
are uniformly distributed and that they track the flow direction accurately. A pulsed laser is set up so 
that the beam forms a sheet of light which passes through the combustor test section perpendicular to 
the direction of flow. The laser light sheet is assumed to be of uniform intensity and thickness. Finally a 
high resolution digital camera is used to capture an image from the laser light scattered by the seed par- 
ticles, resulting in a fuel plume image. The assumptions of uniform particle distribution, particle flow 
tracking, and laser light sheet uniformity are reasonable assumptions which are discussed in more detail 
in the above mentioned references. 


Image Processing 


Fuel plume imaging is a non-intrusive procedure. As a result, the camera view angle is such that the 
raw image will need perspective and optical corrections. As a basis for performing these corrections, a 
grid block with uniform squares marked out on the surface is typically used. This grid block can be 
placed at the same location and orientation as the laser light sheet, and then imaged by the same camera 
which is used to obtain the fuel plume images. By doing this, known image locations from the grid 
block image can be used to properly alter the perspective and optical distortion of the fuel plume image. 

After correcting for perspective and optical distortion, the region of interest must then be extracted 
from the fuel plume image. Ideally the region of interest would encompass the entire combustor cross 
section shape. In reality, laser sheet wall and window interference and reflection patterns can render 
sections of the image unusable. Also, silica particles can agglomerate on the walls and windows creat- 
ing increased or additional interference patterns. Even so, it is possible at times to obtain fuel plume 
images which are clean enough to use the entire combustor cross section as the region of interest. Typi- 
cally an area immediately surrounding a single fuel plume or group of plumes is used as the region of 
interest. 

When the region of interest has been extracted from the perspective corrected fuel plume image, a 
determination must be made of the intensity range which constitutes presence of zero fuel. This is 
accomplished by looking at the intensity values in locations clearly outside of the fuel plume structure 
and identifying the maximum intensity from the ambient background. Once this maximum ambient 
intensity has been identified, all image intensity values in the extracted region of interest are modified 
so that values ranging from zero to the ambient background maximum are reset to a value of zero. All 
other values are proportionally stretched to fill in the range of 1=1 to the maximum image intensity. 
Typical fuel plume images are obtained as 8-bit monochrome images for which the maximum intensity 
range is 1=0 to 1=255. A raw fuel plume image will usually not span the entire range of possible intensi- 
ties, so this image intensity modification will actually increase the apparent dynamic range of the image, 
but this will not affect the correctness or accuracy of the calculated mixing efficiency. A processed fuel 
plume image from which a mixing efficiency will be calculated must have the image intensity range 
modified so that an intensity of zero corresponds to zero fuel. 
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Each individual fuel plume image will have unique subtleties. In some images there may be small 
areas of interference possibly caused by particle agglomeration on the optical window. If these areas are 
small and few, it is possible to utilize either an image filtering process or a photo-retouching process to 
allow the region of interest to include areas which would otherwise be unusable. In some instances, this 
could allow information to be extracted from an image which would have otherwise been entirely unus- 
able. 


Mixing Analysis 


Based on the assumptions outlined above pertaining to uniform particle distribution, particle flow 
tracking, and laser light sheet uniformity, it can be subsequently assumed that fuel plume image inten- 
sity is proportional to seed particle number density, N, when taking particle velocity and image expo- 
sure time into account. 


■image = C 0 NUx where N = ^ pY f (1) 

Cq is a constant along with the image exposure time, T , which represents minimum time duration 
of either the laser pulse, x L ,or camera shutter speed, T c . U represents a known velocity distribution 
across the plane of the fuel plume image. Since a non-constant velocity distribution will affect the 
apparent particle number density, the fuel plume image must be processed to take this into account. 

I = IlM ^ GE so that I = CjN where Cj - C 0 T (2) 

If the velocity distribution is assumed to be constant or is unknown and must be assumed to be con- 
stant, then the following relationships are used. 

I = Ijmage so that I = CjN where Cj = CqUt (3) 

At this point, either equation set (2) or equation set (3) can be substituted into equation set (1) to 
yield the same result. 


I = 


n a 


(4) 


Since both Avogadro’s number, N A , and the fuel molecular weight, Wf , are constant, equation (4) 
can be expressed as follows. 


I = C 2 pY f 


where 



(5) 


Using the average velocity over the region of interest, U , which is constant by definition, yields the 
following. 
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( 6 ) 


I = C 3 pU Yf where C 3 


C_ 2 

u 


Integrating the fuel plume image intensity over the area of the region of interest results in the fol- 
lowing. 

I TO tal = f I dA. = C 3 f pUY f dA = C 3 rii f (7) 

ROI ROI 

Using the definition of the fuel-air mass ratio yields the following relationship. 

m f = f m air = f (pUA ROI - m f ) (8) 

Rearranging equation (8) yields the following. 

(1 + f)m f = f pU A roi and m f = pU A ROI (9) 

Substituting equation set (9) into equation set (7) produces the following. 

^TOTAL = C 3 ril f = C 3rTf pUA ROI = C 3 pU ^f A RO! 

The quantity pU represents the average mass flux over the region of interest, which is constant by 
definition. Next, the average intensity over the region of interest is defined as follows. 

Troi = = C 4 TTf = C A f wtere C 4 = C 3 ^U (11) 

a roi 1 + 1 

Using equation set (1 1) allows the constant C 4 to be redefined using known and obtainable values. 


I roi = C 4 Y, and in general I = C,, Yj- where C 4 = 


l ROI 


1 +f 


( 12 ) 


Next, assume a linear relationship between I ROI and I ST such that at (|) = 1 , I ST = I RO i • 


IrOI _ ° r ^ST _ 


l ROI 


ST 


<l> 


(13) 


The implication inherent in this assumption is that the processed fuel plume image intensities have 
been modified so that 1=0 corresponds to zero fuel. 
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Over the region of interest, for I < I ST , J IdA is proportional to that part of rii f which is fuel lean 
and entirely reactable. For I > I ST in the region of interest, the reactable mass fraction is estimated by 
the following relationship. 


Irx ~ ^4 Y RX — Yf) where Yf — (14) 

Therefore, in the region of interest, for I > I ST , I RX dA is proportional to that part of rfi | which is 
fuel rich and can react only to the extent of available air. Previously [1], this integration has been per- 
formed as a linear approximation between two points, I=I ST and I=Imax- With this formulation, it is 
recommended to perform the integration at each discrete intensity value. One additional consideration 
for this integration is that the upper limit intensity can be no more than a value of I = C 4 . Any fuel 
plume image intensities above this value would numerically supply negative reactable fuel to the inte- 
gration process. This would indicate some non-physical aspect to the processed fuel plume image. 

The definition of mixing efficiency is the following [1,4]. 




JpU Y RX dA 

ril RX 


where Y RX 


f Y f for (j) < 1 
v f ST (l-Y f ) for (j) > 1 


and m RX 


f m f for (|) < 1 
v f ST n \ur for $ > 1 


(15) 


By definition, <|) = f/f ST and f ST = f/([), therefore f ST m air = (f/(j))m air = m f /(|) . Substitut- 
ing this relationship into the appropriate portion of equation set (15) provides the following. 


r rh f for (j) < 1 
m f /(|) for (|) > 1 


(16) 


Using equation set (7) along with the definition of mixing efficiency contained in equation sets (15) 
and (16) allows mixing efficiency to be directly calculated from the information contained in a fuel 
plume image. 




T( | IdA + J ‘RX dA J 

LEAN RICH 

f IdA 

L ROI 

[♦I J IdA + J >RX dA ) 

LEAN RICH 

f IdA 

x . ROI 


for (|) < 1 


for (|) > 1 


(17) 


In arriving at equation (17), note that it is not necessary to be able to obtain a known value for the 
constant C 3 . This is due to the fact that it appears in both the numerator and denominator of equation 
(17). An additional implication of this is that the spanned intensity range of the processed fuel plume 
image does not affect the calculation of a mixing efficiency as long as 1=0 corresponds to zero fuel. 
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Results and Discussion 


The first two fuel plume images were acquired from a Mach 2 scramjet combustor test which repre- 
sented a mid-speed combustor at a flight enthalpy of Mach 5 [5]. These images were both taken at the 
X/H=6 location, where the full test section size is 6H wide by 4H tall. Figures 1 and 2 show the unproc- 
essed fuel plume images, which were both acquired during one test. It should be noted that these images 
do not encompass the entire test section, but do include the entire fuel plume. The result of applying the 
image processing steps to these raw fuel plume images is shown in Figures 3 and 4. These processed 
fuel plume images will be referred to as image 1 and image 2. 




Figure 1. Unprocessed fuel plume image 1 


Figure 2. Unprocessed fuel plume image 2 



Figure 3. Processed fuel plume image 1 Figure 4. Processed fuel plume image 2 


A CFD solution is used to provide velocity field information for use with these fuel plume images 
[5], The velocity field at the X/H=6 location is shown in Figure 5 with the outlined area representing the 
extent of the processed fuel plume images. Figures 3 and 4 contain the result of applying the velocity 
correction to images 1 and 2 respectively. It is found that within the fuel plume, velocity ranges from 
approximately 100 to 1400 m/sec. As a result, it appears that an assumption of constant velocity would 
not work well for mixing analysis of these images. To gain more insight into this matter, it is necessary 
to look at the range of particle travel distances captured in the images. Acquisition of these images 
involved a pulsed laser which formed a sheet approximately 1 mm thick. The limiting factor for the 
image exposure time was the laser pulse duration such that T =T| = 7 nanoseconds. As a result, the par- 
ticle travel distances captured in these fuel plume images range from 0.0007 to 0.0098 laser sheet thick- 
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nesses. Velocity non-uniformity is not an issue in these fuel plume images, which are essentially 
instantaneous images rather than time averaged images. Consequently, even though these fuel plume 
images are taken from a mid-speed combustor with significant velocity non-uniformity, an assumption 
of constant velocity is appropriate. 



Figure 5. Velocity field at X/H=6 with outline of processed fuel plume image area 




Figure 6. Processed fuel plume image 1 with velocity 
correction applied 


Figure 7. Processed fuel plume image 2 with velocity 
correction applied 


Mixing analysis of images 1 and 2 was performed both with and without the constant velocity 
assumption as a basis of comparison between the two methods. For image 1, q M =23.85% assuming 
constant velocity and q M =23.77% with the velocity correction applied. For image 2, q M =23.90% 
assuming constant velocity and q M =23.37% with the velocity correction applied. Applying the con- 
stant velocity assumption, these two images exhibit closer agreement of the resultant mixing efficiency. 

The next set of fuel plume images was obtained from a similar test in the same facility [2], From 
this test, one image taken at X/H=6, image 3, will be analyzed along with an image taken at X/H=2, 
image 4. Computational velocity infomiation is also available for these images. Figures 8-11 contain the 
processed fuel plume images both with and without velocity correction applied. Note that for this set of 
images, the raw fuel plume images did include the entire combustor cross section. And as a result, the 
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processed fuel plume images encompass the entire 6H by 4H test section, a much larger area than that of 
the previous set of processed images. 




Figure 8. Processed fuel plume image 3 



Figure 9. Processed fuel plume image 4 



Figure 10. Processed fuel plume image 3 with velocity Figure 11. Processed fuel plume image 4 with velocity 
correction applied correction applied 


Images 3 and 4 were acquired using a continuous laser which formed a sheet that was 200 p meters 
thick. The limiting factor for the image exposure time for these images was the camera exposure setting 
such that x = T ( . = 0.8 seconds. For image 3, at X/H=6, velocity within the fuel plume ranges from 100 
to 1000 m/sec. In image 4, at X/H=2, velocity within the fuel plume ranges from 1400 to 2200 m/sec. In 
temis of particle travel distances captured in these images, image 3 ranges from 0.4 million to 5.0 mil- 
lion laser sheet thicknesses, while image 4 ranges from 5.6 million to 8.8 million laser sheet thicknesses. 
These images provide time averaged representations of the fuel plume, and as a result, it is desirable to 
apply a velocity correction prior to calculating a mixing efficiency. 

Mixing analysis of images 3 and 4 was also performed with and without the constant velocity 
assumption so that a comparison can be made between the two methods. For image 3, q M =15.02% 
assuming constant velocity and q M =14.02% with the velocity correction applied. For image 4, 
q m =5.63% assuming constant velocity and q M =5.80% with the velocity correction applied. Note that 
at X/H=2 there is less variation of velocity across the fuel plume than at X/H=6, and the resulting mix- 
ing efficiency quantities are closer in agreement at X/H=2 than those at X/H=6. 
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Table 1 provides a summary of results obtained from the mixing analysis of the four images pre- 
sented above. Note that the fuel equivalence ratio is lower for images 3 and 4. As a result, it is not pos- 
sible to make direct comparisons among images 1, 2, and 3 at X/H=6. 


Table 1. Summary of mixing analysis results 


Image 

Type of Image 

X/H 

Velocity 

Correction 

(%) 

1 

Instantaneous 

6 

No 

23.85 

Yes 

23.77 

2 

Instantaneous 

6 

No 

23.90 

Yes 

23.37 

3 

Time Averaged 

6 

No 

15.02 

Yes 

14.02 

4 

Time Averaged 

2 

No 

5.63 

Yes 

5.80 


Conclusions 

When processing fuel plume images to calculate an estimate of mixing efficiency, a known velocity 
field is typically not available. Consequently, the assumption of a constant velocity must then be used. 
For fuel plume images obtained from tests of a mid-speed scramjet combustor, this constant velocity 
assumption appears to be a potential source of great error due to the non-uniform velocity fields present 
in mid-speed scramjet combustors. However, in looking closely at the mixing analysis formulation, it is 
based on integrations of the general fomt j" IdA . Also, as part of the image processing, areas with no 
fuel present will have an intensity value of zero. As a result, an integration of image intensity over area 
will accumulate nothing in areas of zero fuel. Effectively, this means that an assumption of constant 
velocity is really only an assumption of constant velocity within the fuel plume. Non-uniformity of 
velocity outside the fuel plume does not enter into the calculation of mixing efficiency. 

If the velocity field is known when processing fuel plume images, mixing efficiency can be calcu- 
lated both with and without the constant velocity assumption to verify the resulting difference in these 
values. A typical expectation is that the difference between these two values would be greater at loca- 
tions closer to the fuel injectors. Certainly this would seem logical since the non-uniformity of velocity 
is greater across planes that are closer to the fuel injection location. In reality, it does not matter how 
unifomi or non-unifomr the velocity is across the entire plane. Locations closer to fuel injection provide 
greater velocity uniformity across the fuel plume. For mid-speed or high-speed scramjet combustor 
tests, fuel plume image mixing analysis with a constant velocity assumption will be more accurate for 
locations closer to fuel injection. 

It is important to understand whether a fuel plume image represents an instantaneous image or a 
time averaged image of the fuel plume. An instantaneous image does not require a velocity correction to 
be applied since it has effectively been applied already due to the method of image acquisition. On the 
other hand, accuracy of fuel plume image mixing analysis can be improved for time averaged images by 
taking the velocity field into account. This accuracy improvement is greater in situations where the 
velocity non-uniformity across the fuel plume is greater. 
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Summary 


For accurate mixing analysis of fuel plume images obtained from mid-speed scramjet combustor 
tests, two methods can be used to properly take into account the non-constant velocity flowfield. Fuel 
plume images can be acquired as instantaneous images and processed using the constant velocity 
assumption to arrive at an accurate estimate of mixing. And for time averaged fuel plume images, the 
presented mixing analysis formulation for non-constant velocity can be used to obtain an accurate esti- 
mate of mixing if flowfield velocity information is available. 
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